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In nanocomposite magnetic materials the exchange coupling between phases plays a central role 
in the determination of the extrinsic magnetic properties of the material; coercive field, remanence 
magnetization. Exchange coupling is therefore of crucial importance in composite systems made of 
magnetically hard and soft grains or in partially crystallized media including nanosized crystallites 
in a soft matrix. It has been shown also to be a key point in the control of stratified hard / soft media 
coercive field in the research for optimized recording media. A signature of the exchange coupling 
due to the nanostructure is generally obtained on the magnetization curve M{H) with a plateau 
characteristic of the domain wall compression at the hard/soft interface ending at the depinning of 
the wall inside the hard phase. This compression / depinning behavior is clearly evidenced through 
one dimensional description of the interface, which is rigorously possible only in stratified media. 
Starting from a local description of the hard/soft interface in a model for nanocomposite system we 
show that one can extend this kind of behavior for system of hard crystallites embedded in a soft 
matrix. 



I. INTRODUCTION 



A simple model in order to understand the relation between the nanoscale structure and the 
extrinsic magnetization curves in partly crystallized rare earth transition metal alloys (Re-M), can 
be built as an assembly of nanosized crystallites embedded in an amorphous matrix. Such a model 
is based upon experimental caracterization of the samples obtained from an amorphous precur- 
sior, obtained either by high energy milling or a melt spun technique for instance, followed by an 
appropriate annealing which induces the re-crystallization [JQ]. The crystallites and the matrix 
must be characterized by hard and soft magnetocrystalline anisotropy respectively. One key point 
concerning the extrinsic magnetic properties of the system is then the exchange coupling between 
the crystallites and the matrix or between the crystallites via the matrix. In other words, since one 
deals with nanocomposite including magnetic hard and soft phases, such models can enter in the 
so-called exchange spring magnet (ES) systems [s- 11|. Exchange coupling between different phases 



of nanostructred composite materials have been widely studied after the pioneering work of Aharoni 
[12I mainly through the behavior of the average magnetization curve in terms on the applied field. 
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ISHlTl . On the other hand the nucleation field corresponding to the very beginning of 



the magnetization reversal in the soft phase has been determined semi-analytically in different situ- 
ations from the linearization of the micromagnetic energy with respect to the deviation of the local 
magnetization from the direction imposed by the hard phase easy axis [l^ . [igI lisl . [igf . Conversely, 
the local magnetization m(r) in the vicinity of the interface between phases in nanocomposites has 
not been investigated in details but in continuous film geometry (layered systems) 20l-l22l|. where 
furthermore the effect of the hard and soft layers thickness has been investigated through the mag- 
netic phase diagram representing the reversal mode of the demagnetization, namely rigid mode 
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versus ES-couplcd mode, 13- 3] via the analytical determination of the magnetic susceptibility. 
In the works devoted to the magnetization reversal in ES layered systems through the behavior 
of the local magnetization in terms of the distance, z, normal to the interface the solvation of the 
Euler equations relative to the micromagnctic energy minimization is made tractable because of 
the symmetry parallel to the surface according to which the local magnetization depend only on 



20l422l |. This makes the nucleation, domain wall (DW) cornpression and depinning process at 
the hard / soft interface a well established theoretical result 0, 0, 0] which moreover has been 
experimentally observed [s, 23, 24[. In the case of nanocomposites made of true 3D grains either 



simply juxtaposed or separated by grain boundaries, the situation is not so clear and as a general 
rule one refers to the similarity in the AI{H) curve to map the magnetization reversal process to 
the one indeed obtained for the ES layered system. 

In the present work, we deal with a model made of hard inclusions embedded in a soft matrix 
with the purpose to exploit the above mentioned theoretical results on the magnetization reversal 
in magnetically hard/soft composite systems. More precisely our aim is to bridge the gap between 
the stratified media, for which the reversal process is well understood, and our system of magnet- 
ically hard material inclusions embedded in the soft matrix. We investigate the way according to 
which one can map the well-known nucleation, DW compression and depinning process of the mag- 
netization reversal characteristic of the exchange coupled stratified media to the exchange coupled 
inclusion / matrix system. We will use both analytical results for the angle profile at the planar 
interface in a one dimensional approximation and numerical results for true 3-D systems obtained 
from a micromagnetism finite element based code (MAGPAR) |25| . 

In a first step, we compare the 1-D profile to the simulated result for one cubic inclusion embedded 
in a prismatic matrix. The easy axis of the inclusion is parallel to the interface (parallel recording 
media type of geometry). The variation of M{H) in the plateau region characteristic of layered 
ES media is very well reproduced, and the accuracy obtained for the magnetization profile, shows 
that the nucleation, domain wall compression and depinning process is valid for isolated inclusions 
embedded in a very soft matrix. 

In order to investigate the effect of the finite value taken by the edge to edge separation A 
between the crystallites in the actual system, we consider a model made of two cubic inclusions 
in the soft matrix. The magnetization profile behavior in terms of the applied field shows that 
the magnetization reversal follows the nucleation, domain wall compression and depinning process 
when A is larger than the exchange length l^x of the soft material whereas when A gets smaller 
than lf.x the magnetization reverses in a coherent way in the whole system which corresponds to 
the strong coupling regime, or the rigid composite magnet of Refs. 0, [l^l, where the system 
behaves as a single phase. 

Then, we consider the effect of the misallignment between the easy axes of the neighboring inclu- 
sions. In this case, the domain wall which was of Bloch type for the parallel oriented interface, takes 
clearly a Neel character, with nevertheless at the qualitative level, the same type of magnetization 
reversal processes. 
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II. ONE DIMENSIONAL MODEL 



In this section, wc first rccah the basic equations concerning the one dimensional description for 
the magnetization profile accross the hard soft interface. We consider the planar sharp interface 
between a uniaxial hard material and an ultra soft material. The interface, normal to the z axis, is 
defined by z = zq and the easy axis of the hard material, this latter being located at z > zqj coincides 
with the X axis. The micromagnetic characteristics of the media will be denoted Ks.ji, As.ji, Js,h 
where the subscript s{h) refers to the soft (hard) phase and K, A, J denote the magnetocrystalline 
anisotropy constant, the exchange stiffness and /io times the saturation magnetization respectively. 
Here we shall consider the limiting case of an ultra soft phase, with Kg = 0. Then the properties 
of the one dimensional model are totally determined by the reduced parameters sa = Ag/Ah, 
e,7 = Js/ Jh: the reduced distance z* = z/5b; 5b being the Bloch DW thickness in the hard phase. 



6b = HyJ Ah/Kh and the reduced field h — ~Ha/ Hk where Ha and Hk = "^Kh/ Jh are the applied 
and the hard-phase anisotropy field respectively. The minus sign in the definition of h is only for 
convenience in order to deal with /i > in the second quadrant Ha < 0. We consider the case of 
an external field Hq in the x direction. Starting from the energy functional where we explicitly 
assume the magnetization profile to lie in the plane (i, y) and to depend only on z 



E[rh{z)] = S j 




dmx{z)\ f dmy{zy 



dz J \ dz J 
X{z) = XgQ{zo - z) + XhOiz - zo); X = K,AotJ. 



+ K{z){l - m,(z)2) - J{z)m,{z)Ha 



dz 
(1) 



where 0(z) is the Heavyside step function. The minimization of the functional E[rh{z)\ with 
respect to rh{z) leads to the well known Euler equations 
the angular profile (f{z) = {fa{z),x) reads 
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14l I20I - I22I I which, written in terms of 



dz 



±- 



1/2 



with AC7„,,((^(z)) = K{z) (sin2((^(z)) - sin2((^(zf ^^)))) - J{z)Ha (cos(^(z)) - cos((^(z^''"))) 



(2) 



The singularity of the interface leads to the boundary condition at z = zq 



dz 



dip 
dz 



(3) 



In equation ^ the choice for the ± sign depends on the boundary conditions far from the interface. 



represented by the values taken by ipj^ 



s,h 



) where z^ 



denote the location of either 



(sJi) 



±00 we 



the soft or the hard layer center. In the case of the hard / soft interface with z^ 
shall asume in equation ^ that far from the interface, z Zf,, the magnetization 771 is alligned 
to the easy axis and that either Kh^s > or Kh > and Kg = {sin (pi, = 0). Starting from 
an external field applied in the a: > direction, where rh{z) ~ x in the whole system, the field 
is decreased and then the first step in the magnetization reversal is the so-called nucleation field 
h = hnuci corresponding to the reversal in the soft phase, hnuci can be obtained from the expansion 
of E[m{z)] in terms of the small deviations Sm = (m(z) — x) at second order in \Srn\ 18|. In the 
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following wc limit ourselves to a hard layer thick enough for to be considered as infinite. Let 
us first consider the case — > oo. As usual, we consider the situation where the magnetization 
in the soft phase, far from the interface is along the field, namely rh(z ^ •Zb*'') = —x. The profile 
ip{z) can be easyly obtained from the numerical inversion of 



.0 9^ 



dip 



(4) 



which can be integrated, with the result 



In 



1 — tan 



27rVr 



In 



(^/4) \ 



1 + tan((^/4) / 



; for V? > ipn 



A- 



A + cos{(p/2) J 



with A = \ 1- 



Vo 



sin2((^/2) 



for ip < (fo 



(5) 



Here we implicitly asumed that the profile <^(z) is a monotonous function of z and accordingly 
with ipl^'^^ = TT, the soft {z < zq) and the hard (z > zq) phases correspond to tt > ip{z) > ipo and 
(fo > ip(z) > respectively. The value (po of (p aX z = zq is obtained from the boundary condition 

® El 



COs(</Jo) = 



1 



(1 - eAex] 



h{l - eA^j) 



[h\l - CAej f + (1 - eAeK) [(1 - eAex) - 2h[l + e^ej)]] 



1/2- 



(6) 



which is defined up to a critical value which defines the depinning field, i.e. the largest value of 
which can be accomodated by the hard phase, and coincides with the coercive field he since we 
expect the nuclcation field of the soft phase to be smaller, is given by 21 1 



1 - CA^K 

(1 + 



(7) 



is) 

When z^ takes a finite value and then coincides with the half thickness of the soft layer, pt, ~ 
Vbizjj'^^) in equation ^ is considered as a parameter. The profile in equation (0]) is numerically 
integrated and pb is determined from the fulfillment of 



1/2 

'AsiAUsiip^Pb)] dp EE I{pi,) 



with cos{po) = h{l - EAej) + [h'^ {eAej - 1)'^ + 1 - 2h{l - EAe J cos Pb)] 



1/2 



(8) 



from a Ncwton-Raphson procedure by solving : I{pb) — {zjf^ — zq) = 0. 

Let us now consider the case where the easy axis of the hard inclusion is not oriented parallel to 
the interface, {n, z) = Oh 7^ 7r/2, still with (n, y) = -k 12. As a general rule, for (dh — '^/2) not 
too small, the numerical simulations lead to a Ncel domain wall at the interface, and accordingly 
we consider only this situation in the one dimensional model. The magnetization is therefore in 
the {x, z) plane and is totally determined by the angle p{z) = {rh, z). As a consequence of the 
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longitudinal nature of the domain wall at the soft / hard interface the demagnetizing field in the 
soft phase must be introduced. We still do not introduce the demagnetizing field in the hard phase, 
since we aim in fine to model a system where this latter concerns a cubic inclusion characterized 
by Nx = Ny = Nz- The soft phase is enclosed in an ellongated shaped prism, whose long axis, z, 
is normal to the interface. In the following we shall consider the 1-D model in the infinitelly long 
prism geometry with demagnetizing coefficients = Ny = 1/2; A^^ = 0. In the bulk soft phase 
and if z^**^ -> oo, the equilibrium value ip{z^^^) is determined from the minimum of the energy 
density 



£;,, = -J,i7aSin(^<^') + iAiV^sin2((^[^)); with A7V = (TV, - TV,) = i (9) 

Z ^0 I 



with the result 



sin(</j^'^) = - inf (\, ; with : Ha = -Hxh (10) 

In the bulk hard phase, the equilibrium value, Lp"^^ of 'fi{z) is determined in a similar way and can 
expanded in the vicinity of 9h because of the high value of Hk and we get 

1 + asm Uh 

In the present geometry, AUh,si^iz)) are now given by 



h{^{z)) = Ah(-^) (sin2((?„-^(z)) + 2/j(sin((^(z))-sin(^["^))) 



AUsiifiiz)) = A,[—] ■^i^sm\ipiz)) + —sm{ip{z))-sin\^i'')-—sm{^ 



(12) 



where we have introduced the hardness factor k = ^^iqHk I Jh- In a similar way to what have 
been done for the parallel oriented interface we get the profile </?(z) from the solvation of 



Z — Zr, = 



1 /2 

^ (sin^ip) + {2kh/ej) sin((^) - sin\ip['^) ~ i2kh/ej) sin((^[^))) 

(13) 

in the soft phase and 

r*^'^) dip 



Z - Zn ^ - 



,1/2 



J^o (sin^ iip-9)+ 2h{sm{ip) - sin(6'))) 

(14) 

in the hard phase. The value ipo of ip{z = zq) is determined from the boundary condition and 
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results from the numerical solvation of 



sin((po)' + 2h{k/ej) sin((/Jo) - sm{^[''^f - 2h{k/ej) sin((^['^) 



sm{ipo - 9hY ~ sin((^[,''^ - OhY + 2h(s\Q.ipQ - sirn^),"'') 



(15) 



As in the preceding case, when takes a constant value, the value of ip)^ is no more taken from 
the equilibrium condition (|10p but is determined from the fulfillment of z[ip'j^^) = z^"-* with z(<^^'''') 
numerically calculated from equation ()13p . 

Furthermore, we also introduce a way to fit the value of the field, say h^^^^K^ at which the 1-D 
profile is calculated in order to enhance the accuracy with the 3-D simulated profile: instead of 

(s) 

using the actual value of ft,, we fix the value of (^^ to that obtained in the simulation. Hence, 
/j(/»t) satisfies 



^^b^ih) 



(16) 



tfib being the 3-D micromagnetic simulation result for iff,. 

The angular profile across the interface can be used to understand the behavior of the demag- 
netization process in exchange coupled media. This has been done in different situations in the 
literature [a, [igI . [itI . [2^ , especially in the framework of recording media optimisation although the 
magnetization profile accross the interface is generally not explicited. Here, our purpose is to make 
the link with the situation of a lattice of magnetically hard inclusions in a magnetically soft matrix. 
More precisely, we now compare the (f{z) profile calculated from equations (O, (fl^ ITI)) with the 
one extracted from the 3-D micromagnetic simulation in two simple situations. 



III. FEM SIMULATION OF THE GRAIN / MATRIX INTERFACE 

The numerical simulation of the ma gne tization in terms of the applied field Ha is performed 
by the micromagnetic code MAGPAR [25|, based upon a finite elements numerical scheme. The 
embedded hard grain is a cube of edge length a = 2R, with i? = 17 nm being the length scale 
fixed at a convenient value for nanostructured rare earth - transition metal intermetallics [l|-Q. 
We have chosen to keep fixed the magnetic parameters of the hard phase and consequently the 
Bloch domain wall thickness, Sb takes a constant value, 6b = 5nm (see section HVl below) . In the 
model including only one such embedded grain in order to represent the case of isolated crystallites 
in the soft matrix this latter is a parallcpipedic prism of total length = 6R, and lateral width 
Lx = Ly = {2R + S), with 6 ~ 0.40 6b and the embedded grain is located at its center. In the 
model with two crystallites introduced in order to study the influence of A/l^x on the magnetization 
profile, we use a fixed edge to edge distance A = 26b and the ratio A/l^x is varied through the 
value of lex = \/%J^o^7/J^ = 6b / \/ fJt'oiHK / Jh){e.Al ^^j) as a function of ca and ej. The local 
magnetization profile is extracted along a line parallel to the z direction, cutting the grain / matrix 
interface at its center. A particular attention is paid to the quality of the mesh, which must be 
superior to what is necessary for the average magnetization curve, M(Ha) over the whole system. 
Here the quality of the mesh has been controlled through the fulfillment of the unitary condition 
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of m(r) which is exactly satisfied only on the nodes. The typical size of the mesh tetraedra is 
sharply peaked at 0.10 Sb, the maximum edge length for more than half of the tetraedra is less 
than 0.20 Sb leading to a mesh for our (2i? + 6) x {2R + d) x 6R model including about 8.5 10^ 
finite elements. 

Before going further we have to note that because in our system the soft phase surrounds the 
whole hard grain, we have to take into account not only the interface we are interested in, normal 
to z, but also the other sides of the embedded cubic grain. By symmetry, only the top and bottom 
ones are to be considered. This is important especially for the nucleation field determination. To 
estimate the effect of the additional interface, we divide the soft phase domain into the prism 
based on the side normal to z located at z ~ zq namely bounded by {x,y) = ±R, and the top 
and bottom layers {\x\ > R) for which the hard/soft interface is perpendicularly oriented. For 
these additional layers, the demagetizing factors are Ny = — and A^^; ~ 1 and we must add 
to the energy density a term due to the corresponding demagnetizing field and proportional to 
the volumic fraction say a occupied by these layers in the total soft phase domain. Then, in the 
case of the Bloch wall parallel to the interface, the implicit equation from which hnuci is obtained 

Q,H[i3, 

A,, tanh(A,,L,,/2) = e^A, tan{XsL,/2) (17) 



still holds but with modified definitions of the parameters A^ and A,, 



>2 - 



A? = 



OB 
TT 



[1-h] 



N^Js 
t 



(18) 



Since XuLu > 1 for Lh > Sb the value of hnuci is nearly independent of Lh, and therefore, 
wc can estimate the effect of the surrounding layers by solving equation (|17|) for the field h' ~ 
h + aNzJs/ fJ-oHx ■ The nucleation field is thus approximately given by hnuci = h^ — aN^Js/ tJ-oHK 
where ho is the solution of equ. (jl7p with a = 0, and may be negative corresponding to a nucleation 
of the soft phase in the first quadrant. For the same geometrical reason, we do not expect the 
coercive field to be given by the analytical result of the 1-D model (see equ. (O). 

In the following, we mainly focus on the local magnetization profile in term of z for different 
values of the reduced external field, h. 



IV. RESULTS AND DISCUSSION 



We first focus on the magnetization profile obtained at the interface between a cubic grain 
of magnetically hard material and the soft matrix. We consider a magnetically hard material 
characterized by Kj, = 3.05 lO^Jm^'^, ~ 7.7 10^^^ J/m^^ and Jh ~ IT, typical values for the 
Re-Fe compounds leading to Sb = 5nm, the soft matrix being characterized by Kg = and the 
other parameters defined through the values of and ej. Wc first compare the demagnetization 
curves given by the true 3-D model to that one should get from the 1-D profile as calculated 
from equation ([5|). The latter is obtained analytically from the magnetization profile and the 
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geometrical average weighted by the respective volumes of the hard and soft phases in the 3D 
model. The result is displayed in figure ([T]) for 3 characteristic sets of values of (e^,ej). As 
can be seen, the agreement is qualitatively correct, the main discrepancy being the values of the 
nucleation and coercive fields. This stems from shape and finite size effect as as been discussed for 
the nucleation field above. On the other hand, the characteristic plateau in between hnuci and he 
is very well reproduced especially concerning its h— dependence. Then, in figure ^ we compare 
the magnetization profile mz{z) = cos{(p{z)) as calculated in the 1-D description corresponding 

(s) 

to ' = oo to the one obtained from the 3-D simulation. This is done for a set of values 
of h corresponding to characteristic points on the demagnetization plateau {hnuci < h < he). 
We see that on the one hand the compression process of the domain wall in the soft side of the 
interface is clearly evidenced and on the other hand the agreement between the 3-D simulated 
profile and the 1-D anaytical one is quite satisfactory. The only case where the calculated 1-D 
magnetization profile deviates from the simulation result at large values of |z| corresponds to a 
situation {ea ~ 0.75; ej = 0.325 and lex = 2.345 5b) where the plateau mx{z) ~ -1 is not reached 
at the boundary of the system (see figure (Ht)) because the value of the field is too small to compress 
the profile in the limits of the micromagnetic model. We note that the boundary of the system 
corresponds to |z/(5_b| = 6.8. As can be seen in figure ^ the fitting procedure introduced through 
equation (fTB|) definitely solves the problem. We conclude that the 1-D model leads to a rather good 
approximation for the magnetisation curve, and a very good approximation for the angular profile 
at the hard/soft interface. Accordingly the local demagnetization process at the hard grain / soft 
matrix interface follows the one deduced from the 1-D layered model. In particular, we emphasize 
that the strong reduction of the coercive field due to the exchange coupling is reproduced although 
the value of he differs due to the mentioned shape and size effects resulting from the embedding 
geometry. 

Now we focus on the infiuence of the edge to edge distance A between two inclusions, still 
in the parallel orientation, namely n = .t for the two inclusions. The 1-D angular profile is 
calculated by using the numerical determination of ip\^^ as described in section equation ([5]). 
The magnetization profile mx{z) for the value of the field closest to the depinning point of the 3-D 
simulation is displayed in figure ([3]). As ca increases from = 0.162 to 0.75, the ratio {/S./2)/lex of 
the distance between the mid-plane and the interface to the exchange length decreases from 2.115 
to 0.983; hence the fiexibility of the profile mx{z) in the soft phase is sufficient for this latter to 
reach a plateau at 'mx{z) = -1 in between the two hard inclusions only in the first case. Moreover 
when (A/2)/Zex < 1-0 the two inclusions become exchange coupled via the soft phase, and the 
magnetization reversal occurs as a whole with a one phase like behavior, as can be seen in figure 
(g]) where the corresponding demagnetization curve is displayed. This means that in this case, 
a strong coupling reg ime is reached. This is in qualitative agreement with the magnetic phase 



diagram of Ref. [l3| since we clearly evidence the three phases namely the decoupled magnet 



for (A/ 2) /lex > 2, the ES coupled magnet for 2 > {A/2) /lex > 1 and the rigid magnet for 

iA/2)/lex < 1. 

The same conclusions as above hold for the case when one easy axis is not parallel to the interface; 
Here we have chosen as an excmple 6^ = 7r/4; the 3-D simulation has been performed only in the 
two inclusions model with 9h for one of the inclusions. The results for the angular profile 
and the magnetization profiles mx{z) and mz{z) are given in figures ^ and ([7]) respectively. 
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One important difference with the case 6h = is that the magnetization of the hard phase in the 
9h grain presents a reversible variation of rrix before switching as is the case in Stoner Wolfarth 
spherical particles when the easy axis does not coincide with the external field direction. The other 
difference as already mentioned is the fact that the hard / soft interface for h < he is a longitudinal 
(Nccl) domain wall instead of transverse (Bloch) one as it is evidenced from the figure 

Now as we have shown from the local magnetization profile across the interface that the process 
for the demagnetization at the layered system holds at the hard / soft interface of 3D hard inclusions 
in a soft matrix at least for one or two finite sized objects, we relate the demagnetization curve 
of the two inclusions model to an extended one. We consider a lattice including Np = 256 cubic 
inclusions located on the nodes of a simple cubic lattice made of 4 (8 x 8) planes and is an 
extended version of the preceding two grains model. The edge to edge distance between nearest 
neighbors in each plane is A = 26 b as above and the external field is along x. Two different easy 
axes distributions have been used; in the first one all the easy axes arc close to the x -axis, with 
'Ei(\(rii.x)\)/Np = 0.94 and in the second one the axes are randomly distributed on the unit sphere 
{T,i{\{hi.x)\)/Np = 0.5). We have considered the set of parameters as in fig. Qj. Since 2/3~^'^ of the 
interfaces are parallel to both the preferential orientation of the {fii} and the direction of the applied 
field, the comparison with the preceding model for the parallel interface is meaningful. As we can 
see on figure ([5] a), at the qualitative level, the demagnetization curve for the easy axes preferentially 
distributed along x is close the that obtained for the 2- inclusions 'toy' model. The demagnetization 
curve calculated with the randomly distributed easy axes is displayed on figure ([8]b), and leads 
to the same type of conclusion but only at a qualitative level. Given the results we have got on 
the two inclusions 'toy' model, we conclude that the caracteristic plateau in the demagnetization 
curve is indeed the signature of the local magnetization profile compression / depinning process. 
However, considering a mesh of the required quality for extracting the magnetization profile in 
the 256 inclusions model would lead to prohibitively heavy simulations. Indeed, our analysis of 
the local magnetization profile is based upon a continuous function for the later and its behavior 
with respect to the external field. Such a requirement from the finite element simulation results, 
where the profile is obtained from an interpolation between the mesh nodes, can be reached only 
with a very fine mesh. Nevertheless, on the physical point of view, in spite of its continuous 
nature the magnetization profile resulting either from an analytical calculation or a finite elements 
interpolation cannot be interpreted at an infinitely small length scale since the micromagnetic 
formalism remains a continuous medium type of approach valid beyond some characteristic Icnght 
scale, say ^ 1 nni. 

We conclude that as is the case on the small local model, the rigid magnet, exchange spring and 
decoupled magnet regimes are reached on the extended model for roughly A/{2lex) < 1,1 < 
A/{2lex) < 2 and A/(2^ea;) > 2 respectively. In this case, these boundary values can be related to 
the volumic fraction of the embedded cristallites through A/{2lex) = {f v — ^){R/lex) with 
the simple cubic geometry. This extended model provides a link with a more realistic modelling of 
an actual experimental system as was shown in ref. Thus, the results displayed in figure ([HJ 
show that the local demagnetization process explicited at the single intcr-grains soft layer can be 
qualitatively transfercd to a realistic situation. 
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V. CONCLUSION 

In this work, focusing on the well established mccanism of exchange coupling between phases 
in hard / soft composite systems, we made the connection between the demagnetization behavior 
at the layered system where a 1-D description holds and a fully 3-D composite system. The 
connection has been done through the local magnetization profile and we have shown that locally 
at the hard / soft interface this profile can be transfered from the layered to the 3-D systems. While 
it is well known that the distance between hard objects in a soft matrix as measured w.r.t. the 
soft phase exchange length is one of the most relevant parameters driving the coupling, this point 
on the one hand has been made more quantitative and on the hand illustrated by the behavior 
of the simulated local magnetization profile in a small 3-D model. This latter play the role of an 
intermediate between the 1-D description of the interface and the true 3-D model. 
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Figure 1; Demagnetization curve as calculated by using the profile calculated on the 1-D model (lines) 
and compared to the results of the 3-D micromagnetic simulations (symbols). (eA, ej) = (0.75, 0.325) open 
squares and solid line; (0.325, 0.75) open circles and dotted line; (1.00, 1.00) triangles and dashed line. 
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Figure 2: Magnetization profile, mx{z) in terms of z/Sb across the hard/soft interface along the normal 
cutting the center of the cubic inclusion. The hard (soft) phase is located at 2 > (2: < 0). Solid line: 
result of the 1-D model; Symbols: 3-D micromagnetic simulation. Ks = 0; Kh = 3.05 10® Jm'^'^ ; 
a) : = 1, ej = 1 h — 0.0328, squares; 0.0740, circles; 0.1150, upward triangles; 0.1760, downward 
triangles, h) : eA ~ 0.325, e,/ — 0.75 h = 0.0123, squares; 0.1352, circles; 0.2580, triangles, c) : = 
0.75, ej = 0.325 h = 0.0165, squares; 0.0492, circles; 0.1230, upward triangles; 0.2295, downward triangles; 
0.3320, diamonds. Dash-dotted line : result of the 1-D model calculated with z^^'/fe ~ 7.0 and a fitted 
value of the field /i'-''**' = 0.0240, according to equation (fT6|) . 
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Figure 3: Magnetization profile, ma,, (z) between the two cubic inclusions in the soft matrix. A = 2Sb- The 
center of soft layer is located at z/Sb ~ —1.0. Symbols: 3-D micromagnetic simulation; Solid line: profile 
as calculated from the 1-D model. 

a) ej = 0.75 and : = 0.75 and h = 0.1475, triangles; = 0.325 and h = 0.2295, squares; = 0.162 
and h = 0.2622, circles. For ex ~ 0.75, the 1-D model profile is calculated with the fitted value of the field, 
^(/»t) _ 0.1874, according to equation (|16p . In each case, the value of the field is close to the depinning 
field of the 3-D micromagnetic simulation. 

b) ej — 1.50 ; = 0.1623 and h — 0.0697, circes; 0.1885, squares. The larger value of the field is close 
to the depinning field of the 3-D micromagnetic simulation. 



15 




Figure 4: Demagnetization curve of the two cubic inclusions model. e,j = 0.75 ; ea = 0.162, open squares; 
0.325, open circles; 0.75 open triangles. Dashed line : demagnetization curve of the 1-D model for the 
hard/soft interface with e,j — 0.75 and ~ 0.162. 
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Figure 5: Angular profile accross the soft phase between the two hard cubic inclusions located at 
zj&B < —2 and z/5b > respectively with (n, z) = 7r/4 for the latter. Symbols: 3-D micro- 
magnetic simulations for h = 0.0655, triangles; 0.1311, circles and 0.3440, squares. 1-D profile for the 
fitted values of the field /i*^**^ = 0.130, solid lines; 0.2020, dashed lines and 0.3080, dash-dotted lines 



respectively, /i'-'^'*-' is determined using equation H16|) to fit ip{z 



on the simulated result where 



/5b ~ —1 is the location of the soft layer mid plane. 




Figure 6: Same as figure ((5]) for the x— component of the magnetization m{z). 




Figure 7: Same as figure ((SJ for the z— component of the magnetization 
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Figure 8; Demagnetization curve obtained from a micromagnetic simulation on the system including 256 
cubic particles located on the nodes of a simple cubic lattice. 

a) Preferentially oriented easy axes in the direction of the field with Ei(|(ni.a;)|)/A'j, = 0.94 ; ej = 0.75 
and tA ~ 0.162, dashed line; 0.325, solid line and 0.75 dotted line. 

b) Random distribution of easy axes; ej = 0.75 and = 0.162, solid line; 0.75, dashed line. 



